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Abstract. The one variable Krawtchouk polynomials, a special case of the 2^1 function 
did appear in the spectral representation of the transition kernel for a Markov chain studied 
a long time ago by M. Hoare and M. Rahman. A multivariable extension of this Markov 
chain was considered in a later paper by these authors where a certain two variable extension 
of the Fi Appel function shows up in the spectral analysis of the corresponding transition 
kernel. Independently of any probabilistic consideration a certain multivariable version of 
the Gelfand-Aomoto hypergeometric function was considered in papers by H. Mizukawa 
and H. Tanaka. These authors and others such as P. Iliev and P. Tertwilliger treat the two- 
dimensional version of the Hoare-Rahman work from a Lie-thcorctic point of view. P. Iliev 
then treats the general n-dimensional case. All of these authors proved several properties 
of these functions. Here we show that these functions play a crucial role in the spectral 
analysis of the transition kernel that comes from pushing the work of Hoare-Rahman to the 
multivariable case. The methods employed here to prove this as well as several properties 
of these functions are completely different to those used by the authors mentioned above. 

Key words: multivariable Krawtchouk polynomials; Gelfand-Aomoto hypergeometric func- 
tions; cumulative Bernoulli trial; poker dice 
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1 Introduction 

The genesis of this paper goes back to a joint work of Hoare and Rahman [8] in 1983, where 
the idea of the so-called "Cumulative Bernoulli Trials" (CBT) was introduced. The essential 
elements of this probabilistic model are as follows: A player (say, of poker dice) rolls a subset, i, 
of a fixed number, N, of dice for success (say, "aces") with a certain probability a. The player 
is allowed to save his/her k successes, and given a second chance, namely, to mix the i — k 
unsuccessful dice with the previously N — i unrolled ones. The player then rolls the combined 
dice, numbering N — k for success with probability (3. Suppose the number of successes in this 
second try is j — k, which, when combined with the previously earned points, k, gives the total 
number of successes as j. If the number of successes after these two rolls is defined to be the state 
of our system we get a Markov chain by iterating this scheme which took as from state i to state j. 

More explicitly the transition probability matrix of this Markov chain with state space 
0, 1, 2, . . . , iV is given by 



min(i,j) 

K{j,i)= ^ b{k;i;a)b{j -k;N -k;l3), 



(1.1) 
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where 

6(A;;n;p)= (1.2) 

is the binomial distribution. The stationary distribution (po{i) corresponding to this process can 
be defined by 

N 
j=0 

A sufficient condition for some (j)o{i) to satisfy this condition is 

K{i,j)Mj) = K{j,t)Mi), (1-3) 

and also that the summation part of K{j, i) is symmetric in i and j. Use of (|1.2p and (|1.3|) 

gives 

</.o(i) = 6(«;iV;r?), 

where 

(1 - q)^ _ 1 - V _ n-i 
/? "1-/3"''^ ' 

with Di = l + 

^ l—a 

As determined in [8] the eigenvalues Afc of the eigenvalue equation 

N 

'^K{i,j)Tpk{j) = AfcVfe(0, /c = 0,1,2,. .. , 

i=o 

are = a^{l — /S)^, with eigenfunctions 

Mi) = Ki; r]) 2Fi{-i, -k; -N; t?"^) , 

where the hyper geometric function 2-^1 in one variable i is just the Krawtchouk polynomials, 
see |10] . Clearly, by use of the orthogonality property of Krawtchouk polynomials we can write 
down the spectral representation of K{j,i), namely, 

^ /AT\ -1 N/1 0\\ k 



X 2F1 {-i, -k- -N- r?-i) 2F1 {-j, -k- -N- 7]-'). 

It is obvious that the above prototype of dice-tossing and "saving" successes can be extended, 
on one hand, to more practical situations in which Bernoulli trials may be accumulated (for 
example, in 'infection-therapy' models, see [Hl[9]), and on the other hand, to multiple variables 
where one defines various kinds of success (say, aces, kings, queens, ...). Suppose we have 
a process where the total number of dice is A^, of which n subsets ii, . . . , i„ are tossed separately 
for successes of n different kinds with probabilities ai, 02, ■ ■ ■ , an, respectively. Let ki,k2, ■ ■ ■ ,kn 
be the number of successes in each category (that is, ki aces, k2 kings, . . . , etc.). Mix the 
"unsuccessful" ii — ki,i2 — k2, ■ ■ ■ ,in — kn trials with the remaining N — ii — 12 — ■ ■ ■ — in dice. 
The player "saves" the k successes and is allowed to try again with the resulting N — ki — ■ ■ ■ — kn 
dice with probabilities /3i, /32, • • • , /3n of producing ji — ki aces, j2 — k2 kings, . . . , etc. 
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With the definition for the n-fold multinomial distribution 6„ given below we get that the 
transition probability kernel from the state (ii, ^2, • • • , in) to (ji, j2, . . . , jn) is 

K{ji,...,jn;ii,---,in) (1-4) 

n 

= ^ • • • ^^n(jl - fci, . . . ,jn - kn;N - ki- k2 kn;(3l,...,f3n)Y[ KK^ir] OLr) . 

ki k„ r=l 

The stationary distribution in this case, as before, is defined by 

^(j;i)0o(i) = i^(i;j)</'o(j), i = (n,Z2, . . . ,^n), j = (ii,j2, • • • , jn), (1-5) 
which, together with (ll.4p gives, as a sufficient condition 

^'i'---''"^ fc=i \ j=i 

which is the n-fold multinomial distribution. Using ()1.4p and (jl.Sp we find that, in perfect 
analogy to the one-variable case, the r/'s are related to the probability parameters a's and /3's 
in the following way: 



1 - E^fc 

??i = — V2 = ■■■ = — rjn = = 1 ~ ^'^'^^ = ^' (^•^) 



1 - "1 1 - a2 1 - an fc=i . ^-i 
i]i = — r]2 = ■■■ = — = ^ " 

Pi P2 Pn 



1 - E A fc=i 

fc=l 



where 



I?n = l + VT^^, 0<afc, /3fc<l. 

One of the questions we will address in this paper is this: what are the eigenvalues and 
eigenfunctions of i^(i;j)? In other words, we will seek solutions of the eigenvalue problem: 

N 

^(i; j)V'k(j) = AkV'k(i). (1.7) 

jl,...,jn=0 

In the single- variable case the eigenfunctions are simply 0o(i) times the ordinary Krawtchouk 
polynomials, as we mentioned earlier. So it is reasonable to expect that in n dimensions (n > 2) 
the solutions will be an appropriate extension of these polynomials. The question is: which one? 
Even in the n = 2 case the question is not quite as straightforward as it would seem. In fact, 
the same authors, Hoare and Rahman, wrestled with this problem for a number of years until 
they were able to show that the Krawtchouk limit of the 9 — j symbols of quantum angular 
momentum theory in physics provides the answer which, written in slightly more convenient 
notation is 

(2) 

b2{xi,X2;N;r]i,r]2)Fl ' {-mi, -m2; -xi, -X2; -N;t,u,v,w), (1.8) 

(2) 

where {t,u,v,w) satisfy certain relationships with 7]i and r]2, and the represents an iterate 
or a 2-variable extension of the familiar Fi Appell function: 
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that is, 

k j k I ''■J-'^-''-\^)i+]+k+l 

subject to conditions for convergence in case they are infinite sums. Unbeknownst to the authors, 
Hoare and Rahman, at the time they pubhshed the paper |9], a general n variable extension of 
the ^ in ([L8]), namely, 

n n 

\[{-mi)n \[{-Xi)^ 

F;"'(-m,-.-iv;u):= E — nM-A^)g." n-S'. (1-^) 

a special case of Gelfand hypergeometric function [2l H], was known to, among others, the 
Japanese authors Mizukawa and Tanaka [15], who used them to prove their orthogonality in 
some special cases. 

Later, Mizukawa [13] gave a complete orthogonality proof of p.9|) using Gelfand pairs, 
then [15] used character algebras and closely following this proof came Iliev and Terwilliger's 
proofs, first for n = 2 [TT], then for general n in [12], in which they used tools from Lie algebra 
theory. In these proofs the authors found it convenient to use 1 — Uij instead of Uij as parameters 
in ()1.9|) . and also to use 

n n n 

r]Q = l-'^r]i, mo = N-'^mi, xo = N-'^Xi. (1.10) 

i=l 1=1 i=l 

Our second objective in this paper is to give an alternate proof of orthogonality by using the 
more elementary method of hypergeometric functions and their various transformation proper- 
ties and integral representation, and in doing so we find no special advantage of using 1 — Uij 
instead of Uij, or of using (jl.lOp . Also, we believe that elementary and perhaps a bit cumber- 
some as it may be, our method yields a byproduct that throws some light on the underlying 
geometrical structure of these polynomials, which the other authors may have overlooked. 

In Section [2] we will list the transformation and integral representation formulas for the 
multivariable hypergeometric functions that we shall use throughout this paper. Section [3] will 
be devoted to obtaining the necessary conditions of orthogonality, while in Sections H] and [5] 
we shall deal with the sufficient conditions that will simultaneously establish the orthogonality 
relation 

^6„(x; N- r7)Ff "X; -N; u)F^") (-m'; -x; -N; u) = 5^,„,V^n(m; iV; r?), (1.11) 

{x} 

where rj are the parameters for the dual orthogonality. The relationship between the n's and 
the 77's and 77's will be given in latter sections, specially Section [6l 

In Section[6]we shall examine the geometrical implications of the relationships between the u's 
that result from the necessary and sufficient conditions, while the last section will be aimed at 
proving that the n-dimensional extension of the function in (II. Sp are precisely the eigenfunctions 
of the kernel i^(x;y). In Section [7] we get the expressions for the eigenvalues and eigenfunctions 
of the transition probability kernel governing the evolution of our Markov chain. 

Before closing this section we must mention that the polynomials in (|1.9p were, in fact, 
introduced into the statistical literature by R.C. Griffiths [5] as early as 1971 which he defined 
as coefficients in an expansion of their generating function. However, Mizukawa and Tanaka jl5j 
seem to have been the first to give the explicit expression in (jl.9p . 
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Few people would disagree with the importance of looking at certain mathematical objects 
from different points of view. We feel that this is certainly valid in the case of the present prob- 
lem: character algebras, Gelfand-Aomoto functions, Lie algebras and the much older methods 
of hypergeometric functions including their series as well as their integral representations have 
a useful role to play. Such a wealth of approaches may be important if one tries to obtain matrix 
valued versions of these probabilistic models in the spirit of [BJ [7j . For a very rich and recent 
extension of the scalar valued solution of the hypergeometric equation, see [16]. It is worth 
noticing that in this case there is yet no extension of the Euler integral representation formula 
for the 2-Fi function, and that algebraic methods such as those coming from Lie algebras or group 
representation theory, which have played such in important role in [IH \T7[ [T3t [T^ [T5] have not 
made a mark in this matrix valued extension yet. Similar consideration would be relevant if one 
were to consider a matrix valued extension of, for instance, the work in [3]- 

2 Transformation formulas and integral representations 
(a) Transformation formulas 

The point of this section is to establish (j2.ip . (j2.2p and (j2.3p . 

For references to some of the classical identities in this section the reader can consult for 
instance [1]. 

For n = 0,1,2,..., 

2Fi{—n, a; c; x) = 2-?^i(— a;l + a — c — n;l — x). 




By one iteration 



Fi{-n;a,b; c;x,y) 



{c — a — b) 



n 



n; a,b;l + a + b — c — n;l — x,l — y). 



(c). 



By multiple iteration 




JTT-l, . . . , TTT'ni Xl, . . . , Xyj, 



TV; nil, • • • ,Uln,U2l, ■ ■ ■ ,U2n,Unl, 





i=l 




F 



,(n) 





i-N) 



1 



n 



1=1 



the second identity follows provided the x's are nonnegative integers with ^Xi < N. 
The second transformation for the 2-^"! function is: 
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whose first iteration gives 

Fi (a; b, b'; c; x, y) = (1 - y)-''Fi (a-b,c-b- b'; 

V l-y l-yj 

= (1 - (a-c-b- b', 6'; 

\ 1 — X 1 — X 

and the n-th iteration gives 

= (1 - (1 - n2„)-"^ •••(!- «nn)-"" 

X ai,...,a„;/3i,...,/3„_i,7-/3i /3„;7; 



«11 - ^^In ^^12 - Uin -Uln Uin - Unn -U' 



nn 



1 - Uln - Uln 1 - Uln 1 " Unn 1 - U, 

which is vahd for ah a's, /3's and 7, provided the series remain convergent. 

(b) Integral representations 

For < Re 6 < Re c, 

2Fi(a, 6; c; x) = ^...^^ ^'"^1 " - Cxy^d^, 

T{b)T{c - b) Jo 

Fi{a;b, c; d;x,y) 



(2.2) 








r(6)r(c)r(fi-6-c) ^ 

o<5i+6<i 

provided Re(6, c, d, d — b — c) > 0. 
In n dimensions this extends to 

F^")(ai,...,a„;/3i,...,/3 

n ) T) ^11 ) • • • 5 ^In ) • • • ) ''^nl ) • • • ) ''^rm ) 



r(7) 



n \ rt 



1 fl 



7- 



j=i \ k=i 

provided 

< Re ai, . . . ,an,7 - 



i=i J 



This is the formula that we shah find most useful throughout the paper. 
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3 Necessary conditions of orthogonality 

The point of this section is to show that (j3.4p is a necessary condition to insure the desired 
orthogonality (j3.2|) . 

Denoting the F^"^^ polynomials in (jl.Op by Pni(x) for abbreviation, we find that 

n \ 

> (3-1) 

i=i / 

which is just the generating function of Pm(x). In a sense this represents the opposite point 
of view of Griffiths [5] where he defined the polynomials as the coefficients of the generating 
function. 

Since one of our aims is to prove the orthogonality relation 

C' :=5^^'n(x;A^;r7)P^(x)P„,(x) = if m / m', (3.2) 

{x} 

it must follow, as a necessary condition, that 

J^6„(x;7V;t7)P^(x) =0, m/(0,...,0), 

{x} 



which, by (13. ip . implies that 

n \ 

i=i / 

and therefore 

n 

'^r]jUij = l, i = l,2, ...,n. (3.3) 

i=i 

If we denote the parameters of the dual orthogonality by = {rji,r]2, ■ ■ ■ ,r]n)^ then, using 
5„(m; N;ri)P^{^) = 0, x (0, 0, . . . , 0), 

{m} 

we get as a necessary condition 

n 

^f]jUji = l, i = l,...,n. (3.4) 

4 Sufficient conditions of orthogonality 

In this section we show that the relation (|4.6p among the parameters is sufficient for orthogo- 
nality. 

Instead of the sum in p.2p let us consider, for the time being, the following sum 

^6„(x;iV;77)Pm(x)F|"^(ai,a2, • • • -mi, -m„;7;uii, . . . ,uin, . . . ,'u„i, . . . ,'u„n), 

{x} 
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in which we assume that while the m's are nonnegative integers the a's are not, and nor is 7 
a nonpositive integer, and that 



< Re Qi, ■ ■ ■ ,an,7 - 



ai 



1=1 



Using (j2.3p we then find that the above sum equals 

n 



r(7) 



n \ " 

r ( 7 - Ea* nr(ai; 



1 i-\ 



JO 

n 

o<EeKi 



i=l 



n er"'^?^ E r,)p^(x) n h - E 

{x} j=l \ 



i=l 



However, 



TV 



j;6„(x;iV;77)J] 1- J^e.n. 

{x} j=l V i=l 



n u ~ E ^^^^^ 



k=l 



i=l 



(4.1) 



But, by (j3.3p . ^ ^j'^ij = Ij so ^CiVjUij = ^ ^j, and hence the sum inside the integral in ()4.ip 

j=l i,j i=l 

becomes 



1-E& E- 



E ■'•ij 



-E^j 



i=l 



i,j k=l 



1 - E ^^^if" 



E 



4=1 



r=\ 



where 



Ar = E'?«'"rs - 1. r = 1,2,... ,n, 



s=l 
n 



^r,s — ^ ^ 'Hj^Tj'^sj 1 — ^s,r5 r, 5 — 1, 2, . . . , 



We now expand the n-fold product on the right-hand side of (|4.2p to get 



En 



rrii 



. 1 \kil,ki2,---,kin—lj . 

,}l=l ^ ' ' ' ^ 1=1 



n-1 



J = l ^fcll+fc21^fci2+fc21 . . . ^'=1,1 



n-1 

mi- E ki,j+{k2,i+k„^i) 

w /lfc22+fc32 /lfc23 + fc42 Akn,n-l+kn-l,n-l ^ 3 = 1 

^ ^2,3 ^2,4 ■ ■ ■ ^n,n-l si 



(4.2) 



(4.3) 
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n — 1 71—1 
Yl fc2,j+(fcll+^32H \-kn2) rrin- Yl kn.j-\-ikl,n-l+k2,n-l-\ \-kn-l,n-l) 

x^, . (4.4) 

Substitution of ()4.2p and (j4.4p inside the integrand of (j4.ip and computing the integral, with 
the Qj's replaced by — m^'s and 7 by — A^, gives the value of the integral as 

r— 2 J — 1 r — 3 J — 1 

X • • • X (-m^) „_i (m'l H h - A)Ar_,mi m„ 

r = l j = l 

n-1 

X^fcn+/=21...^fc^n_-l+^n-l.n-l-Q^^^ . = 1 ^(_iV)^. (4.5) 

r=l 

Let m = ej, m' = e^- , where ej = (0, . . . , 1, 0, . . . , 0), = (0, . . . , 0, . . . , 1, . . . , 0), i.e., e's are 
unit vectors with the i-th component being 1 for and j'-th position for e^-, otherwise all 
components 0. 

In particular, for m = ei, and m' = e„, 

n-1 

1- E fci: 



1 / 1 \ 1- '^'ij , 

N^{^[ku,...,h,^J^' ^^^^^^ 

{'=i,il 

X(0) n-1 (0)fcn(0)fcl,2---(-lk.n-l =-%^- 

1- E fcij ^ 
j=i 

By similar arguments it follows that 
Ai,j = 0, j, 
which, by (|4.3p . means that 

n 

^r/jUr-jtisj = 1, r, s = 1,2,... ,n, r / s. (4.6) 



So the only terms that survive in (j4.5p give 

+ E rn'j „ n {-< 

k=l / N-V m, k=l 



kJ^k 



n 



I™' = TT A""^ (4 7) 



The first factor in the numerator on the right-hand side of (j4.7p is zero unless A — ^ > 

A:=l 

n 

A - E^T-fc) i-e-, 
fc=i 

n n 

^m, >j;m;., (4.8) 

i=i j=i 

while the n remaining factors imply that they are zero unless 
fn[ > rrii, i = 1,2, . . . ,n. 
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(|4.7p and ()4.8p imply that I™' = unless m' = m. So we have 



ruk n 



fc=i 

and that (j4.6p is a sufficient condition for orthogonality. Clearly, it can be written in a form 
closer to the desired form (|l.lip . namely 



TV ^ ^m,m' • 



( ) 



(4.9) 



In the following section we complete this part of the work by expressing Y[ terms of 

the parameters of dual orthogonality, i.e., ?7j's. 



5 Reduction of dlTQl ) to ( fTTB 

As a typical A let us consider 

n 

Ai = Y,risul-l by (113]). 



From (|3.3p we have 

n 

^r/j^xij = 1, 
while (14.61) gives 



(5.1) 



rjjUijUrj = 1, r = 2, . . . , n. 
From dO and (1521) it follows that 



'71 



^12^13 ■ ■ ■ ^Ir 

A 



11... 1 

1 U22 ■■■ U2n 

1 Un2 ■■■ Unn 
1 

Ull 



Mil . . .Ui^r-lUl,T+l ■■■Uin 

A 



knl 



where 



k=l 



1 1 

U21 U22 

Unl Un2 



1 



1 1 

U2,r~l 1 
Un,r—1 1 



1 

U2n 



k=l 



(5.2) 



r = 2,...,n. 



D = Det (uij), 
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which follows from (13.31). So 



,s=l 





1 


1 . 


1 






1 


1 . 


1 


Uu 


1 


U22 


U2n 


+ Ul2 


U21 


1 . 


U2n 




1 


Un2 
1 .. 


1 


1 






1 . 


Unn 


+ Uln 


U2l . . 


• U2,n- 


-1 1 




-Dm 










Unl . . 




-1 1 











+ 



(5.3) 



It is clear that expansion by the first rows of the n determinants on the right produces n 
terms: 



D + uu 



1... 1 

1 U22 • • • k2n 
1 Un2 ■■■ Unn 



+ U12 



1 1 ... 1 

U21 1 U23 ... U2n 



-\ '^Uin 



1 

U21 

Unl 



1 

Un,n—1 1 
Un,n—1 1 



whence, we derive by careful rearrangement of the determinants the following expression: 



D 



Uu Ui2 ... Uir, 
1 1 ... 1 



""31 U32 ■ 
Unl Un2 ■ 



UZr, 



nil Ui2 



Ulr 



Un-1,1 Un-1,2 ■ ■ ■ Un-l; 
1 1 ... 1 



From ()5.3p and ()5.4p we finally obtain the desired expression 



(5.4) 



E2 1 s=l 
rjsUis - 1 = = , 

s=l '1 

and, in a similar way, we derive the general form: 



E2 1 s=l 

Vr 



r = 1, 2, . . . , n. 



6 Relations among the parameters 

The previous sections were devoted to proving the orthogonality of the n^-parameter polynomials 
Pm{x) defined in (|1.9p . with respect to the multinomial distribution 6n(x; N;r]),r] = {r]i,...,r]n) 
(or 6n(m; N;r]),r]= (r/i, . . . , r/„), with E = E ^j)- Clearly, the n relations in (13. 3p completely 
define the r/j's in terms of the u's, as they do the r/j's, by (j3.4p . If we think of the r/'s as given 
parameters then the polynomial parameters must satisfy in addition — (2n — 1) = (n — 1)^ 
relations among themselves. Where do they come from and what do they mean geometrically 
is the question we shall examine in this section. 
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Note that the sufficient conditions of orthogonaUty (j4.6p give (2) further relations between 
the T/'s and n's. In fact, combining (|3.3|) and (|4.6|) we get 

'^r]jUrj{l - Usj) = 0, r,s = 1,2, ... ,n, r ^ s. (6.1) 

i=i 

Clearly, the positivity of all the probability parameters, i.e., r/'s, require that the determinant 
Det(itrj(l — Usj)) = 0. It is possible to analyze these determinants and obtain the missing 
(n — 1)^ relations. We shall, however, take a different route. 

Our proof of orthogonality is based almost entirely on the integral representation (12. 3p and 
the simple multinomial summation formula. The proof is direct and elementary. However, if 
we had instead used the transformations (j2.ip and (12. 2p we could have reduced the problem to 
a (n — l)-variable case, with the polynomials having (n — 1)^ parameters, instead of n'^. There 
are, of course, n ways to make this reduction, depending on which of the n variables we integrate 
first. If we do the Xn summation first, then by n long and tedious set of computations we can 
find that the reduces to an expression proportional to the (n — l)-fold sum 

/■X1J-X2 /-X„-l 




X\, . . . , Xji—l, ^ 1 



UnnUll , UnnUU 



. , UnlUln UinUn2 

1 ,1 ,...,1 ,...,1 

UlnUn2 UinUn,n—l 'U\nUn,n—l 

(1 - n„„)(l - uii) 




-m„_i; -xi, ... , -Xn-i; -^mj; 1 



.^-^ (1 - n„i)(l - 

^ {1 - Unn){l - Uln) ^ (1 - U^n) (1 - ) ^ (1 - ttnn) (1 - ^21 ) \ 

(1 - ?Xln)(l - ?in2) ' " ' ' (1 - ni„)(l (1 - ti2n) (1 " 'Wn2) ' " ' 

where 

r]iUnl{l - Unl) . mUn2il-Un2) . 7/„_i'Un,n-l (1 

U — r; T' ^2 — T, C,n 



If we were to have orthogonality of the n-variable case then surely a reduction to a lower 
dimension will retain the same property. However, the appearance of the two F^^ functions 
above doesn't suggest their parameters are the same. But the point is that they are, not 
identically, but consistent with the necessary conditions of orthogonality (13. 3p . To illustrate 



this point let us take the product of ^'s and the first n — 1 parameters of the first -F^^" ^\ i.e. 
compute 

Unl-UinJ V UinUn2j \ UinUn,n-l 

n—1 

^ _ 7?iUii(l - Unl) _ 772^12(1 - Un2) _ _ - Un,n-l) 

n-1 

n-1 J2 Vj'^lj ~ '^nj) n-1 , ^ n-1 

2^Ci — r"2^'^* + Tvz. — T"2^^* + ^ by dsu). 

Vn'^lnyl- l^nn) Vnl^nnyl- '^nn) 
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On the other hand, 

(1 - n„„)(l - un) 



Ci 1 



(1 - «nl)(l - Uln) 



(l-nnn)(l--»12) ^ , 

+ <,2(1 - — — -) H 



(1 - ui„)(l - n.„2)' 



+ Cn-1 1 



(1 'Unn)(l ^n,n2) 
(1 - Uin){l - Un,n2) 



n-1 

EC. 

i=l 



??l^»l(l - Uii) 
'UnUnniX ^In) 



??2Mn2(l - ^12) 



(1 - Ui^n-l) 
VnUnni^ - Ui^n) 



n-1 



Similarly the equalities of other parameters are also established. The (n — 1)^ relations can be 
expressed in the following compact form: 



(6.2) 



where 



(one must, of course, take it for granted that Uij ^ 0, for any i, j). 

The geometrical implication of these relations seems to suggest that the n-variable Kraw- 
tchouk polynomials (jl.9p live on an (n — l)^-dimensional submanifold, defined by (|6.2p . of the 
n^-dimensional space of (uij). 



7 Eigenvalues and eigenf unctions of K{[;]) 

For the sake of notational consistency and convenience we shall replace K{i;j) by K{x;y), so 
that by use of (II. 4|) and (II. 7p . our eigenvalue problem can be stated as follows: find ^m(x) such 
that 



6„(x; N; (3) bniy; N- r?) J] " E 

{y} k=l {r} 

^1 / ^, \ 



ai 



a. 



PniOn - 1) 



(-^l)ri ■ ■ ■ i-Xn)r„i-yi)ri " " " (-yn)r 

ipmiy) = Xmbn{x;N;r])ip^{x.). 



(7.1) 



Mm - 1) 

We will show that 

V'm(x) = Pm(x) 

for some choices of the parameters Ujj's. To this end what we will do is compute the sum above 
with ^m(x) replaced by Pm{^)- 

Using the integral representation (j2.3p we find that the sum that needs to be computed is 



n "1 

{y} k=l 3=1 I 



j=i I V i=i 



i=l j=l \ i=l 



i=l 



A'- En 



(7.2) 
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Since by (|1.6p 



a, 



1 - E n 

1 - E A ^=1 

«=1 



we have 



a,- 



1=1 

n 



Therefore the right-hand side of (|7.2p simphfies to 



(-iV) 



riH \-rn 



i = l 



i=l 



yl-Eft , 
\ i=i / 



1 - E ^ft 



1=1 



1 - E * J 
j=i 



n 



TV 



(-AT) „ 

1 1-Eft 
\ i=i 

n f / n 

xll 

i=l I \ 1=1 



1-Eft 
1 - E * 



1 - E 

J 1 - E A i 



So the sum over the rj's is 
N 

i=l 




n / \ n 
Xk 



1=1 



j=i / j=i \ J i=i J j=i 



n 1 ^ " E ^^''ii ~ " E ^^^^ 
j=i I 



j=i 
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1 - ^ ^iii 



i=l 



with 



■jUij 



Since 



1=1 



n 



the right-hand side of (|7.ip . via (j2.3p . becomes 



5„(x;iV;77) 



1 rl 



JO 

o<E 5«<i 



i=l 



j=i [ i=i J \ j=i 



n(-mi)„+i n(-^i)", 

j=l Y.kiji=l Z^,'=J» 



n%!(-iV) n „+l 

i=ij=i 



n n+1 



X (Xi H h Xn - A^) 



J]4-"+i J] J]{a^.n,, + (1 - a.Oa;,}'^-- 

i = l «=1 J=lj = l 



N 



i=l E'=iji=l E 



ajUjj + (1 - gjQcjt - 



1 - Wi 



So Pm(x) is an eigenfunction of A'(x;y), provided Uij exists for all z,j = 1,2,., 
that 



(7.3) 



. , n, such 



aAuij - uoi) = (1 - iOi)uij, 



I.e., 



Q;j(l - (jJi) - aj{l - Uij) = (1 - uji)uij, 
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^zii-M = l_^^ = l_Y^ - (7.4) 

So 

aj{l-Uij) _ Uij{l-aj) 



PkUik = 1 - /3jUij — - — = -PjUij + 

OLj Uij 

l3juf, - (1 - Qj + ajl3j)u 



aj Uij 

On the other hand, from ()7.4p it follows that 



i = 1, 2, . . . , n. 



a,(l-Ui,) Qfc(l-'Uifc) ak{l-aj)ui 

~ ^ Uik — 



Oj-Uij Oik-Uik aj{l - Uk) - {aj - ak)uij 

So the equation to solve for Uij, j = 1,2, . . . ,n, for each i, is 

(1 - Oj) _ a , ,^ „, \ "fc/^fc 



— /?,• + (1 - a,) — — , 

Uij-aj aj{l - ak) - [aj - ak)Uij 

since we must have Ujj ^ 0, which has n roots. Detailed analysis of these roots is not of 
immediate interest to us. 
For n = 2, 

1-ai 02/32(1 -ai) ai/3i(l - a2)+ "2/32(1 - a2)/3i(ai- a2)Mii 

~ Pi + ~ 



nil— ai ai(l — 02) — (ai - 02)^11 ai(l — q;2) — (ai — 02)^11 

which leads to 

/3i(ai - a2)iuii - ai)^ - (1 - ai)(ai - 02 + ai/3i + a2/32)('Wii - "i) + ai(l - ai)^ = 0, 

which is exactly the same as formula (5.7) in [9]. 

Finally, it is clear from (j7.3p that the eigenvalues are 



Xm = n(i 



i=l 
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